SECOND ORDER DIFFERENCE EQUATIONS AND DISCRETE 
ORTHOGONAL POLYNOMIALS OF TWO VARIABLES 



YUAN XU 

Abstract. The second order partial difference equation of two variables 

Vu := A m (ie)A 1 Viu + A li2 0r)AiV 2 u + A 2 ,i(x)A 2 Vim + A 2t2 {x) A 2 V 2 u 

+ B 1 (x)A 1 u + B 2 (x)A 2 u = Xu, 

is studied to determine when it has orthogonal polynomials as solutions. We 
derive conditions on T> so that a weight function W exists for which WDu is 
self-adjoint and the difference equation has polynomial solutions which are or- 
thogonal with respect to W . The solutions are essentially the classical discrete 
orthogonal polynomials of two variables. 



1. Introduction 

Let A and V denote the forward and the backward difference operators, respec- 
tively, 

A/(t) = /(t + l)-/(t) and V/(t) = /(t)-/(t-l). 
The second order difference equation in one variable takes the form 

(1.1) a(t)AVu + b(t)Au = Xu, Ael. 

In the case that a and b are polynomials, it is known that some of the difference 
equations can have orthogonal polynomials as solutions. Such difference equations 
have been characterized completely and the orthogonal polynomials so derived are 
the classical discrete orthogonal polynomials (Hahn, Meixner, Krawtchouk and 
Charlier polynomials), see pQ for a survey in this direction and for references. 

The purpose of the present paper is to consider the second order difference 
equations in two variables. Throughout this paper we use the notation 

ei = (l,0) and e 2 = (0,l). 

For x = (xi,X2) € K 2 we define the forward and the backward partial difference 
operators by 

(1.2) A iU (x) = f(x + <a) - f(x) and V iU (x) = f(x) - f(x - *), » = 1,2, 
respectively. We consider the second order partial difference equation 

(1.3) Vu := Ai,i(a;)AiViu + A 1 , 2 (x)A 1 V 2 u + A 2 ,i(x)A 2 Viu + A 2 ,2(x)A 2 V 2 u 

+ B 1 (x)A 1 u + B 2 {x)A 2 u = Xu, 
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where Ai t j and Bi are polynomials and A is a real number. The goal is to charac- 
terize those equations that have orthogonal polynomials as solutions. 

The orthogonal polynomials that satisfy l|1.3|) will be discrete orthogonal poly- 
nomials of two variables. A sequence of polynomials, {p a }, of d variables is called 
a sequence of discrete orthogonal polynomials if there exist a set of isolated points 
V C R d and a weight function W defined on V such that p a are orthogonal with 
respect to the bilinear form 

xev 

Since we will discuss orthogonal polynomials associated with the difference equation 
(|1.3|) . the set V will be a lattice set, that is, a subset of Z x Z. We should mention 
that a simple change of variables x — ► hx shows that we can replace the difference 
operators Aj by the divided difference 

[f(xi + h,x 2 ) - f(xi,x 2 )]/h and [f{x 1 ,x 2 + k)-f(xi,x 2 )]/k. 

In this case the second difference AjVj becomes second order divided difference; 
for example, A1V1 and A1V2 become 

[f{xi + h,x 2 ) - 2f(xi,x 2 ) + f(xi - h,x 2 )]/h 2 and 

[f{xi + h,x 2 ) - f(xi + h, x 2 - k) - f(xi,x 2 ) + f(xx,x 2 - k)]/hk, 

respectively. In other words, the difference equations l|1.3|) is the model for the 
uniform grid. 

After the dilation indicated above, letting (h,k) —* (0,0), the divided difference 
operators become the differential operators; in particular, Aj becomes d/di and 
A^Vj becomes d 2 /didj. This way the difference equation becomes a differential 
equation. In this connection, it should be pointed out that the classification of 
the second order partial differential equations that have orthogonal polynomial as 
solutions was carried out in j^j, see 00 El El El for further work in this direction. 

Comparing to the theory in one variable, the structure of discrete orthogonal 
polynomials in several variables is much more complicated. Some basic results are 
obtained in JS]; the relevant ones will be recalled in the following section. 

In the next section we will also discuss for what A{ _j and Bi the equation (|1.3() 
has polynomial solutions. Such equations are called admissible. In Section 3 we 
define and discuss the self-adjoint form of the second order difference equations 
and give necessary and sufficient conditions for the existence of a function W so 
that WD is self-adjoint. We also define a notion of W being consistent with the 
difference equation, and show that consistence and self-adjointness will imply that 
the difference equation has polynomial solutions that are orthogonal with respect 
to W. In Section 4, we identify difference equations that have discrete orthogonal 
polynomials as solutions. Up to the mild restriction that we put on the difference 
equations, we believe that these are the only such difference equations. The poly- 
nomials are the classical discrete orthogonal polynomials of two variables (see, for 
example, jUElE]) on the set V = {x : x\ > 0, x 2 > 0, x\ + x 2 < N} and on the 
set V = {x : < xi < M, < x 2 < N}, where both N and M can be infinity. 

2. Preliminary and Admissible equations 

2.1. Discrete orthogonal polynomials. Let V be a subset of R d . Orthogo- 
nal polynomials on V depend on the structure of the polynomial ideal I(V) := 
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{p G R[ir] : p(x) = 0,Va: G V} that has V as its variety. The discrete orthogo- 
nal polynomials on V can only consist of polynomials that do not belong to I(V). 
Such polynomial subspaces have been identified and the basic structure of dis- 
crete orthogonal polynomials of several variables, including three-term relation and 
Favard's theorem, has been studied in |15|. 

Let M.[V] denote the subspace of polynomials. This is the space that orthogonal 
polynomials belong to. Then M.[V] = R[xi, . . . , Xd}/I (V). There is a lattice set 
A = A(V) such that every polynomial P G B>[V] can be written as 

P(x) = c a x a mod I{V), c a G M. 

aeA 

One particular result in shows that the set A satisfies the following property 

a G A implies a — f3 G A, whenever a — (3 G Nq and (3 G Nq. 

For two variables, d = 2, this means that A must be of a stair shape; more precisely, 
there is a sequence of positive integers rij, which satisfies n rn < n m _i < . . . < no 
(some of the rii can be positive infinity and so can m), such that 

(2.1) A= {(k,l) : < I < m,0 < k < ni}. 

Hence, any polynomial P G M.[V] can be written as 

m ni 

We will assume that A contains {(i, j) : < i+j < 2} and at least one of the points 

(2.1) and (1,2). 

Given A as in l|2.1|l . the highest degree of monomials appeared in R[V] is denoted 
by n, which could be positive infinity. It follows that 

n := max{n; + 1 : < I < m}. 

For A = A(V), denote Ah(V) — ■ i + j = k} for each k < n. Let denote 

the distinct elements in Afc(V^), which is the number of monomials of degree exactly 
k in R[V}. For each j satisfying < j < n, it can be verified that 

m 

(2.2) r k = k + \-Y,{k-l-n l ) + , 

;=o 

where (a)+ = if a < and (a)+ = a if a > 0. 

Let W be a real function on V and W(x) ^0on7 for any x G V. Assume that 

K x i\\W(x)\ < oo for all G Nq 

in the case V is an infinite set. Let II 2 — H[xi,X2] and let II 2 denote the subspace 
of polynomials of degree at most n. Define the bilinear form (•, •} on II 2 x IT 2 by 

(f,9) = C(fg), where £(/):= £ f(x)W(x). 

If (/, g) = 0, we say that / and g are orthogonal to each other with respect to W on 
the discrete set V. A polynomial of degree k is said to be an orthogonal polynomial 
with respect to W if it is orthogonal to all polynomials of degree lower than k 
in M.[V]. For a given weight function W, let Vk denote the space of orthogonal 
polynomials of degree exactly k on V . The dimension of Vk is in 1)2. 2J1 . Only in 
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the case of V — {x : x\ > 0, x 2 > 0, < x\ + X2 < N} this dimension is ru = k+1, 
which is equal to the dimension for the space of continuous orthogonal polynomials 
(see0). 

Since we are interested in discrete orthogonal polynomials that are solutions of 
the difference equation l|1.3|l . the set V is a lattice set. If u satisfies the difference 
equation then the function u* = u(x\ — k,x 2 — I) will satisfy a similar difference 
equation. Hence We can assume that (0, 0) belong to V. For a given G ZxZ, 
let Bij denote the square that has (0,0), (i, 0), (0, j) and as corners. In order 
that the differences make sense, we further assume that V satisfies the following 
property: 

(2.3) (i, j) 6 V implies (k,l)&V, whenever (k, I) € Bij. 

As pointed out in JS], in some cases, the index set A can be the same as the set 
V. This happens, in particular, in the following two important cases, 
Type A. V = {{x 1 ,x 2 ) : < x x < M, < x 2 < N}, 
Type B. V = {{x 1 ,x 2 ) ■ Xi > 0, x 2 > 0, < x x + x 2 < N}, 

that are of particular interest to us. We note that both N and M can be infinity. 

2.2. Admissible difference equations. Since the forward and the backward dif- 
ference equations are related by 

(2.4) A^VjAj+Vj, 

a second order partial difference operator can be expressed in different forms. We 
shall fix one particular form and call it the standard form. 

Definition 2.1. Let Ai_j and Bi, i,j = 1,2, be polynomials of two variables. We 
call the operator T> defined by 

(2.5) V :=A M (x)AiVi + Ai, 2 (x)AiV 2 + ,4 2 ,i(z)A 2 Vi + A 2:2 (a;)A 2 V 2 

+ B 1 {x)A 1 + B 2 {x)A 2 
a standard second order partial difference operator. 

An alternative way of writing a second order partial difference operator is 
V :=1i,i(x)AiVi + !i, 2 (a;)AiV 2 + A 2 ,i(x)A 2 Vi + l 2 , 2 (x)A 2 V 2 
+ B 1 (x)V 1 + B 2 (x)V 2 , 
which can be easily converted to the standard form by using (|2.4Jl . 

Definition 2.2. The equation T>u — Xu is called admissible on V if for any k 6 No, 

there is a number Afc such that the equation T>u = XkU has r^ linearly independent 
polynomial solutions and it has no non-trivial solutions in the set of polynomials of 
degree small than k. 

The following proposition shows that we can assume that the linearly inde- 
pendent polynomial solutions are given by those polynomials whose highest term 
contains a single monomial. 

Proposition 2.3. The equation T>u = Xu is admissible in R[V] if and only if for 
every integer k that satisfies < k < n, there exists a number such that the 
equation T>u = X^u has r^ linearly independent solutions in the form of 

(2.6) Pkt{x)=xU l 2+Rk,i(x), (M)eA fc , 
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where Rk,i S n^, and the equation has no non-trivial solutions in K[A] of degree 
less than k. 

Proof. That the equation is admissible implies that, for each k, there exist solutions 
Qij E K[V], degQij — k. Introducing the notation 

x\ = {x\x{:{i,j)eK, i + j<k} 

and regard it as a column vector, we can write 

Qij(x) = afjxX + R hj , degRij < k 

where aj ■ is a row vector of real numbers. We show that {aj \jx\ ■ (i,j) E A} 
is linearly independent. Indeed, suppose otherwise, then there exist Cjj such that 
Y^u,j)eA Ci J a T,j x A = 0, which implies that 

y c i,jQi,j = : ^ i- 

(i,i)GA (j,i)eA 

Since ^ satisfies Vu = XkU, we conclude that 

VRk = V ^ CijQij = Xk ^ CijQij = XkRk 
(i.i)eA (ij')eA 

so that Rk = by the definition of admissibility. Consequently, since the set of 
{Qi,j} is linearly independent, Cjj = for (i,j) E A. Thus, {afjX^ : (i, j) E A^} is 
linearly independent and the matrix {&J j){i.j)£A k is invertible. Let 

(-^'^(ij'JeA* : = [( a I,j)(i,j)eA k ] (Q»,j)(i,j)eA fc ; 

then Pj ; j are polynomials of the form Pij(x) — x\x 2 +pij where degp^.j < i+j = k. 
This proves one direction of the proposition. The other direction is trivial, since 
the set {Pi,j} in l|2.6|) is evidently linear independent. □ 

Proposition 2.4. For the equation T>u = Xu to be admissible it is necessary that 
Aij are polynomials of degree at most 2 and Bi are polynomials of degree at most 
1. ' 

Proof. Suppose that T>u = Xu is admissible. Then there is a system of polynomials, 
P{i,j} as in (|2.6() . which are solutions of the equation. Let A®j and B® denote the 
part of Aij and Bi that has degree larger than 2, respectively. We need to prove 
that A®j(x) = and Bf(x) = 0. Substituting P^j into Vu — Xu and examing the 
highest terms of the resulting equation, we obtain 

A\±(x)A x Vxx\x{ + A° 12 (x)A 1 V 2 x\x 3 2 + A° 21 (x)A 2 V 1 x\x 3 2 

+ A° 22 (x)A 2 \7 2 x\x 2 + B°(x)A lX \x 2 + B° 2 {x)A 2 x[x{ = 

for (i,j) E A. Let = (1,0) and = (0,1); then AiV jx\x J 2 — and we 

conclude that Bf(x) = B 2 (x)° = 0. Next we let = (2, 0); then only AiVitffa^ 
is nonzero and, consequently, A\ ^(x) = 0. Similarly, we deduce A 22 (x) = by 
setting (i,j) — (0, 2). Furthermore, the choice (i,j) = (1, 1) leads to 

<2(z)+4i(*) = - 

Finally, let (i, j) = (1,2) (or (2,1)); we conclude that 

(2xi + l)A h2 {x) + (2xi - 1)^,1 (at) = 0. 
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Together, the last two displayed equations imply A® 2 (x) — A\ 1 (x) =0. □ 

Since Ay are polynomials of second degree and Bi are polynomials of the first 
degree, we can assume that they have the form. 

A f \ 2,0 2 , 1.1 , 0,2 2 , 1.0 i 0,1 i 0,0 

A ij i x ) = a % \j x i + ij x ^ + a i,j x 2 + «ij x i + a i,j X? + a i,j ; 

Bi(x) = bjxi + b 2 x 2 +b° l , i,j = 1,2, 

Theorem 2.5. Assume that the quadratic parts of A\ )2 and A 2t i are equal; that is, 
Ax t 2 — A2.1 G 11^. Then the equation T)u = Xu is admissible if and only if 

(2.7) Ai t j(x) = axiXj + a^jXx + bijx 2 + Ci.j and Bi = bxi + di 
for i,j = 1,2 with ci,i + ci,a + c 2 ,i + c 2 , 2 = and 

(2.8) A fc = k(ka -a + b), < k < n, 

where a and b are real numbers such that ap + b ^ for all nonnegative integer p 
satisfying < p < n. 

Proof. Again consider the highest terms of the polynomials in T>P kt i — \Pk,l- Since 
Aix{ = (xi + l) 3 - x\ = jx\~ + . . . and Vix\ = x{ — (xi - l) j = jx\~ + . . ., we 
obtain 

A 1 s(x)k(k - l)xi~ 2 x l 2 + Ai :2 (x)klx\~ 1 x 1 ^ 1 + A 2A (x)klx k l ~ 1 x\~ x 

+ 1(1 - l)A 2 , 2 (x)x k x l 2 2 + B l (x)kx\^ 1 x\ + B^x^x^f 1 = X ]+k x k x l 2 . 

Using the formula of Ai j and Bi, we can then derive equations on the coefficients 
of Aij and Bi by comparing the coefficients of x\x\. 

The coefficients of x\~ 2 x l 2 2 and x k+2 x l 2 ~ 2 lead to equations 

k(k - 1)0^1 = and 1(1 - l)a 2 2 '° = 

which imply that \ ~ a 2 ' 2 — 0. Using the assumption that a k ' l 2 — a 2 'i for 
k + I = 2, the coefficients for x\~ 1 x l 2 +1 , x\x l 2 and x\ +1 x 2 _1 lead to the equations 

k(k - l)a}{\ + 2Ma°{ 2 2 + kb\ = 0, 

(2.9) k(k - 1)0,1'° + 2kla{' 1 2 + 1(1 - l)a 2 \ + kb\ + lb\ = X k+U 

2klal'° 2 + 1(1 - l)a\^ + lb\ = 

for (k, I) e A. In particular, if (k, I) — (1, 0) and (k, 1) — (0, 1) then the first and the 
third equations of l|2.9|l show that b\ = and b\ = 0, respectively. Furthermore, 
k = 1 in the first equation shows that a{ 2 — and, consequently, a x \ = 0. 
Similarly, I = 1 in the third equation shows that a 2 ' 2 — and a 22 = 0. Setting 
k = or I = in the second equation of i|2.9|l leads to the equations 

1(1 — l)<z 2 '2 + lb 2 — Xi and k(k — l)ct^'° + kb\ — X k . 

In particular, k = I = 1 gives b\ = b\ and, consequently, a 2 {\ = a 22 . Furthermore, 
letting (k,l) = (1,1) and (k,l) = (2,0) in the second equation of l|2.9(l gives the 
equation 

2a\ 2 + 2b = X 2 and 2al'" + 2b = X 2 
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where we write b = b\ = b\. Consequently, we conclude that = a\'\. Writing 
a = a^'J = a\'2 = 0,221 we have proved that Aij and Bi are in the forms specified. 
Moreover, the second equation of (|2.9|l with I = or with k — shows that 

fc(fc — l)a + kb = X k and 1(1 — l)a + lb = Xi. 

This is (|2.8fl and that A; 7^ shows that ap + b ^ for all nonnegative integer p 
that satisfies < p < n — 1 . 

On the other hand, assume that Aij and Bi are given as in lj2.7[l and Xk is given 
as in (|2.8|) . Let = tit + 1) . . . (t + k — 1) be the Pochhammer symbol. Define 

m k (t) = {t) k /k\ and m k j{x) = m fe (xi)m ; (x 2 ). 

For each pair of nonnegative integers (r, s) satisfying r + s = n, substituting the 
polynomial 

P r , s (x) = m r ,s{x) + ^2 fkjm k ,i(x) 

k+Kn 

into the difference equation T>u = X n u leads to a linear system of equations for the 
coefficients fij. Indeed, the definition of m k j implies that 

fc-i 

Am k (x) = m k -i(x + 1) = ^m^^x) and Vm k {x) — m k -i{x), 

j=o 

from which it follows that 

fe-2 fe-1 

AiVim fe j = y^TOj,;, and AiV 2 m fe j = ^m^y-i, 

and similar formula for ViA2m& i i and A2V2?n/ £i i. Furthermore, there is also 
tmk{t) = (k + l)mfc + i(i;) — km k {t), from which follows formulas such as 

x\&.xV\mk,i(x) = k(k - l)m k ,i - (k - l) 2 m k -i,i 
xiX2AiVim k j(x) = klm k j(x) - k(l - l)m,k,i-i(x). 

Using these formulas, a tedious computation shows that the coefficients of m r ^ s in 
the expression T>P r s — X n P r , s are linear in f k j with r + s < k + l < m. Since the set 
of m k j is clearly linearly independent, these coefficients must be zero, which gives 
a linear system of equations in variables f k ,i- 

The coefficient of m r , s turns out to be c\ t i + c\p, + + 02,2, which is zero 
by assumption. Clearly, this also shows that this condition is necessary. The 
coefficients of m r —i tS and m r>s -i lead to the following two equations 

(r - l)(-a(r - 1) + 01,1) + s(-a(r - 1) + a 2 ,i) + s&i,2 + di + K-i,sU-i,a = 0, 
r(-a(s - 1) + ai, 2 ) + (s - - 1) + 02,2) + rb 2 ,i + d 2 + X r , s -if r , s -i = 0, 

from which / r -i, s and /r,s-i are determined uniquely. The coefficient of m^; for 
fc + I < r + s is a linear combination of / p g in which p + q > k + l and is the 
only element with p + q = k + I. The coefficient of f k ,i is 

X k ,i - X r:S = (k + l-r-s)(b + a(-l + k + l + r + s)) 

which is nonzero by assumption. Hence, starting from / r , s — 1 and f r ,s—i and works 
our way down, we conclude that f k j is uniquely determined. □ 
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If 11 satisfies an admissible difference equation T>u, then under a change of vari- 
ables x i— > x + h, where h 6 Z 2 , the function u(ai) = u(.t + ft,) will satisfies a 
similar admissible difference equation. We consider these equations equivalent. In 
particular, we can choose h such that B±(x) — bx± and B2(x) = bx 2 . 

We note that an affine transform xi = pi,ix[ + pi, 2 x 2 , x 2 = Pi.\x\ + P2,2x' 2 , 
however, will alter the difference equations significantly. 

The characterization of partial differential equations in two variables that have 
orthogonal polynomials as solutions is given up to affine transformations; in other 
words, if u(x) satisfies a differential equation, then the function u*(x') := u(x) 
under the affine transformation also satisfies a similar differential equation and the 
two equations are considered to be the same type. This allows us to use affine 
transformation to reduce each type of equations to its simplest form and proceeds 
from there. For the difference equations, however, the affine transformation can no 
longer be used as a tool to reduce the type of equations. 

3. Self-adjoint difference equation and the weight function 

3.1. Self-adjoint difference equations. For the difference equation in one vari- 
able, the self-adjoint form is 

A(a(x)Vu(x)) - Xu(x) = 0. 

For difference equations in two variables, we introduce the following definition: 

Definition 3.1. Let a^j be continuous junctions of two variables. The operator 

£u = Vi [cti,i(x)Aim + a\ >2 (x)A 2 v\ + V 2 [a2,i(x)A 1 u + ct 2i2 (iz;)A 2 m] 

is called a self-adjoint second order difference operator. The equation £u — Xu is 
called a selj- adjoint difference equation. 

Proposition 3.2. Let T> be the difference operator in (|2.5() . There is a function 
W such that W{x)T> is self-adjoint if and only ifW satisfies 

WB X = Ai(WA 1A ) + A 2 (WA h2 ) 

(3 ' X) WB 2 =A 1 {WA 2 ,i)+A 2 (WA 2 , 2 ). 

Proof. First we rewriting the self-adjoint form of £u in the standard form of the 
difference operator T>u. For this purpose the following formulas are useful: 

(3.2) A i (fg) = fA i g + g(- + e i )A i f and S7 i (fg) = fV i g + g(--e i )V i f 
where i = 1, 2 and 

(3.3) V i u{- + e j ) = AjViU + ViU = AjS/iU + AiU- AiViU, i,j = 1,2; 
all can be easily verified. Using these formulas we can write 

£ u =cti i i(- — ei)AiViw + <ti,2(- — e 2 )AiV2M 

+ c 2 a(- - ei)A 2 Viu + er 2:2 (- - e 2 )A 2 V 2 u 
+ (Vi<ti,i + V 2 (Ji i2 )Aiu + (Vitrei + V 2 ct 2 , 2 )A 2 u. 
Comparing with the standard form T>u in (|2.5|l and solving for er^ leads to 
= o"i,i(- - ei), A 1>2 = cri )2 (- - e 2 ), Si = Vcr^i + V 2 eri,2, 
^2,1 = 02,i(- - ei), A 2 , 2 = o- 2 . 2 {- - e 2 ), B 2 = Vi(Ti, 2 + V 2 a 2j2 . 



DIFFERENCE EQUATIONS AND ORTHOGONAL POLYNOMIALS 



9 



Consequently, it follows that 

= + e i)i l<*,j<2. 
Using the fact that Vjit(- + ej) = Ajtt, we then obtain 

Bx = AiAi,i + A 2 Ai !2 and 5 2 = AiA 2ll + A 2 A 2 , 2 - 

That is, the operator T> is self-adjoint only if the above relations hold. Evidently, if 
B\ and -B 2 are given by the above relation, then T>u can be written in the self-adjoint 
form. Hence, T> is self adjoint if and only if these relations hold. Consequently, WT> 
is self-adjoint if and only if (|3.1|l holds. □ 

Using (|3.2|) , the system of equations (|3.1[) can be written as 

Ai A (x + ei)W(x + ex) + Ax, 2 {x + e 2 )W{x + e 2 ) = n(x) 

[ ' >A ' A 2 ,x{x + ex)W(x + ex) + A 2 , 2 (x + e 2 )W(x + e 2 ) = r 2 (x) 

where are polynomials defined by 

n = Ax t x + A lt 2 + Bx and r 2 = A 2 j + A 2i2 + B 2 . 

If the determinant of this system of equations is nonzero, then we can solve for 
ViW and V 2 W. Let us define 



ct{x) 
(3x{x) 



=A lt x(x + ex)A 2 , 2 (x + e 2 ) - A ia (x + e 2 )A %l {x + e x ), 
=ri(x)A 2:2 (x + e 2 ) - r 2 (x)A lt2 (x + e 2 ), 
=r 2 (x)Ai^x(x + ex) - rx(x)A 2tl (x + ei). 



Theorem 3.3. Assume that a(x) ^ for x G V. Then there is a W such that 
WT>u is self-adjoint if and only if 

(3.5) /3x(x)a(x + e 2 )f3 2 (x + ex) = f3 2 (x)a(x + ex)(3x(x + e 2 ), x E V. 

Proof. Since a(x) ^ means the determinant of the system of equations (|3.4() is 
nonzero, we can solve (|3.4|l to get 

( 36 ) W(x + ex) _ f3x(x) ^ W(x + e 2 ) _ (3 2 {x) 



W[x) a(x) W{x) a[x) ' 

These equations imply 

rir / . . \ W / . ,P2(x + ex) J3x{x) f3 2 (x + ex) 

W{x + ex+ e 2 ) = W(x + ex) , — r- = W(x — -— — — - — f, 

a(x + ex) a(x) a(x + ex) 

itt, \ Txr/ ,fli(x + e 2 ) TT _, . (i 2 (x) 3x(x + e 2 ) 

W(x + ex+ e 2 ) = W(x + e 2 ) ' ' ' = W(x) " ' ' ' ' ' 

a(x + e 2 ) a(x) a(x + e 2 ) 

Consequently, we conclude that 

(3x(x) p 2 (x + ex) = (3 2 (x) (3x(x + e 2 ) 
a(x) a(x + ex) a{x) a{x + e 2 ) 

from which the equation l|3.5|) follows. 

On the other hand, assume that l|3.5|l holds. We need to show that the system of 
equations l|3.4(l has a solution. We denote 7j :— f3j/a. The condition (|3.5(l becomes 

/ 3 7) li(x-e 2 ) = 72(3; -ei) 

71 (x) 72 (x) 
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Recall that the set V satisfies 1)2. 3|) and (0, 0) 6 V. Clearly, multiplying 71 or 72 
by a nonzero constant will not change the equation, so that we can assume that 
7i(0) = 1 and 72(0) = 1. For x\ > and 22 > 0, we define a weight function W by 

W(x) =71 (2 -e 1 )7i(x-2ei)... 71 (0,2 2 ) 72 (0,ar 2 - 1)72(0, x 2 - 2) ... 72(0). 

Using the equation 1)3. 7J) . we then have 

ty(a; + ei) = 71(g) 72(1,22 - 1) 72(1,1)72(1,0) 
W(x) 71(0, 12) 72(0, x 2 -l) "'72(0,1)72(0,0) 
71 (2) 71(0,^2) 7i(0,2)7i(0,l) 
71(0,22)71(0,2:2 - 1) "'7i(0,l)7i(0,0) 
=7i 

Furthermore, using the equation (|3 - 1 Of) again, 

W(x + eg) = 7i(gi - 1, x 2 + 1) 71 (an - 2,2 2 + 1) 71(0, 2 2 + 1) 72(0, 2 2 ) 
W(x) j 1 (x 1 -l,x 2 ) 71(21 -2,2 2 ) "' 71(0,^2) 72(0,0) 

72(2) 72(2 -ei) 72(1,22) 72(0,22) 
72(2 - ei) 72(2 - 2ei) ' ' ' 72(0,22) 72(0,0) 
=72(2). 

Consequently, the function W defined above is a solution for the equations l|3.4[l . □ 

3.2. Another notion of self-adjointness. We should like to point out that, just 
as in the definition of the standard difference equation, we can reverse the role of the 
forward and the backward differences and give another notion of self-adjointness. 

Let ciij be continuous functions of two variables. We could call a second order 
difference operator £ self-adjoint if 

(3.8) £u — Ai [5i,i(x)Viu + 01,2(2;) Vait] + A 2 [o 2 ,i(2)Vitx + a2 : 2(2)V 2 u] . 

For this notion of the self-adjointness, we can also determine the conditions under 
which WV is self-adjoint. Let us define 

5(2) :=(A lt i(x - ei) + B 1 (x))(A 2i2 (x - e 2 ) + B 2 {x)) - A 2 ,i(2 - e 2 )A lf2 (x - ei), 
n(x) :=A 1A (x) + A 2 ,i(x), and r 2 (x) := ^1,2(2) + ^2,2(2). 

Following the proof that leads to the result for V, we can prove the following: 

Proposition 3.4. Let T> be the difference operator in (12.51) . There is a function 
W such that W(x)T> is self-adjoint in the sense of <|3.8I) if and only ifW satisfies 

WB 1 = V 1 (WA ltl )(- + e 1 )+V 2 (WA lr2 )(- + e 1 ) 
WB 2 = Vi(WA 2tl )(- + e 2 )+V 2 (WA 2 , 2 )(- + e 2 ). 

Furthermore, assume that a(x) ^ 0; then there is a function W such that WT>u is 
self-adjoint if and only if 

2(2 - ei)r 2 (2)Fi(2 - e 2 ) = a(x - e 2 )rx(x)r 2 (x - ei). 



DIFFERENCE EQUATIONS AND ORTHOGONAL POLYNOMIALS 



11 



3.3. Weight function and orthogonality. Given a weight function W defined 
on a lattice set V, we define 

= £ mw( X ). 

xeV 

For the difference operator given in D, we define 

(3.10) Ax = A 1;1 Vi + A ia V 2 + Si and .A 2 = A 2 ,i Vi + A 2a V 2 + B 2 . 
From the definition of T> we have 

(3.11) Vu = A1A1M + A 2 A 2 u. 

Let dV denote the boundary of the set V; that is, dV is a subset of V such that 
if x G <9V then at least one of the elements x ± +ej, i = 1,2, does not belong to 
V. In order that the difference equation makes sense on V, it needs to satisfy some 
condition on dV . For example, we can impose the conditions Q3.12[l below on the 
coefficients Aij. 

Proposition 3.5. If WD is self-adjoint and Aij satisfy the conditions 
Ai i(x) = A 2 i(x) = 0, x G dV and x + e x G V 

(3 12) 

Al,2(a;) = ^2,2(2:) = 0, leWandi + ejeV, 

then, for any polynomial g, 

C{Aig) = and C(A 2 g) = 0. 

Proof. Let ftV = {x G <9F : x + e x G V"} and d 2 V = {x £ dV : x + e 2 £ V}. 
Changing summation index and using the fact that W(x + e\) = if x £ <9V^ \ 9i^, 
we obtain 



£(Ai,iVis) = J! M / (^)Ai,i(^)Vi5W 



= - ^ Ai(WM M )(x) 5 (x) + £ T^(x)Ai 4 (x) 5 (x). 

A similar equation holds for C(Ai^ 2 V 2 g). Hence, by the assumption on Aij, we 
conclude that 



C(A ltl V l9 ) = -J2 MWA 1A )(x)g(x), 
xev 

£{A lt2 V 2 g) =-J2 A 2 (WA h2 )(x)g{x). 



xev 

Since WD is self-adjoint if and only if l|3.1|l holds, by Proposition 13.21 it follows 
that 

C{Aig) =£(Ai,iVi 5 ) + £(A 1>2 Vig) + C(B ig ) 

= -J2 [Ai(Wi4i,i)(a:) + A 2 (WA 1<2 ){x) - W{x)B 1 {x)\ g[x) = 0. 
xev 

The proof that C(A 2 g) — is similar. □ 
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Definition 3.6. The weight function W is consistent with the difference operator 
V if for all g 6 II 2 , 

(1) C(Aig) = and C(A 2 g) = 0; 

(2) C{A h2 g{--e 2 ))=L{A 2 , 1 g{--e 1 )). 

In particular, setting g{x) = cc, shows that if W is consistent with T> then 

C{B X ) = and £{B 2 ) = 0. 

By Proposition ^ .51 if WD is self-adjoint and the coefficients of V satishes H3.12I) . 
then W is consistent with V if the condition (2) in the above definition holds. We 
need the additional condition for the following result: 

Theorem 3.7. If W is consistent with the difference operator T> and the polyno- 
mials u and v satisfy T>u = Xu and T>v = \xv, then 

(X-fi)C(uv) = 0. 

Proof. By the definition of Ai and the product formula for A.,- , 

Axtfg) = A 1A [fV ig + g{- - ei )Vi/] + A 1<2 [fW 2 g + g{- - e 2 )V 2 /] + BJg 
= fAig + <?(• - e 1 )A ljl V 1 f + <?(■ - e 2 )A 1 , 2 V 2 /. 

Applying C to both sides and using the fact that W is consistent with V, we 
conclude that 

C(fA l9 ) = -£(<?(■ - ei)Ai,iVi/ + <?(• - e 2 )A 1)2 V 2 /). 

Similarly, working with A 2 , we have 

C(fA 2 g) = -C(g(- - ei )ViA 2)2 Vi/ + <?(• - e 2 )^ 2 , 2 V 2 /). 

From these two identities and the fact that Aju(- — ej) = Vjit, it follows that 

CivAiAw) = -£(VimAi,iViu) - £(Aiu(- - e 2 )A h2 W 2 v), 
C(vA 2 A 2 u) = -£(A 2 u(- - ei)A 2 ,iViu) - C(W 2 uA 2t2 W 2 v). 

Hence, by . 1 1|> and the fact that Aju(- — ej) = VjU, it follows that 

(A - fi)C(uv) =C(vAiA-lu) + jC(vA 2 A 2 u) - C(uAiA-lv) - £(uA 2 A 2 v) 
= - £(Aiu(- - e 2 )A h2 V 2 v) - C(A 2 u(- - ei)A 2 ,iViu) 

+ C(A lV (- - e 2 )A ia V 2 u) + C{A 2 v{- - ei)A 2 ,iVitt) 
= - C{A x , 2 {A x uA 2 v){- - e 2 )) - £(A 2A (A 2 uA lV )(- - e x )) 

+ C(A h2 (A lV A 2 u)(- - e 2 )) + £(A 2A {A 2 vAiu){- ~ d)) = 

by the part (2) in the Definition 13.61 of the consistency. □ 
Corollary 3.8. If WD is self-adjoint and the coefficients of T> satisfies l|3.12[l and 

(3.13) C(A h2 g(--e 2 )) = C(A 2A g(-~e 1 )), jetf, 

then the equation Vu = X^u has polynomial solutions that are orthogonal with 
respect to W on V . 
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3.4. Condition for consistency. Let us look at the conditions for the consistency 
closely. Recall that the Type A set V — {(xi, x 2 ) ■ < xi < M, < x 2 < N} and 
the Type B set V = {(xi,x 2 ) ■ £1 > 0, x 2 > 0, < x 1 + x 2 < N}. If V is of Type 
A or Type B, the condition l|3.12[l means 

(3.14) i4i,i(0,X2) = i4a,i(0,aj a ) = 0, and Ai, 2 (xi,0) = A 2>2 (x u 0) = 0. 

Since A\ t \ and A 2 ,i are coefficients of A1V1 and A2V1, both involve Vi, and A 2t \ 
and A 2 .2 are coefficients of A1V2 and A2V2, both involve V2, the conditions in 
(|3.14|l appear to be reasonable assumptions. 

Recall that d x V = {x G dV : x + e x G V} and d 2 V = {x G 8V : x + e 2 G V}. 

Proposition 3.9. Assume that A\ y2 — A 2: \ G If|. If T> is admissible and (|3.12|) 
holds, then d\V is a subset of {x : xi = 0} and d 2 V is a subset of {x : x 2 — 0}. In 
particular, (|3.12(l becomes ()3.14|l . 

Proof. The condition 13.12f) means that the pair Ai ! i(x) and ^2,1(2;), and the pair 
^1,2(2;) and ^2,2(2;), agrees on the set dV\. Recall that both A\ t \ and A 2i i are 
quadratic polynomials and they are of the form l|2.7() . If dV± contains more than 
4 points, then the classical Bezout theorem on the algebraic curves will imply that 
Ai t i and A 2 ^ has a common factor. Since Ai^{x) = ax\ + 7Ti ! i(x) and A 2i i(x) = 
ax\X 2 + 7T2,i(a;), where TTij are polynomials of degree 1, the common factor has to 
be a linear polynomial, which means that both polynomials are product of linear 
factors. Thus, they must be of the form 

Ai,! (x) = ax\ (xi + 61) and A 2>1 (x) = ax\{x 2 + b 2 ). 

Consequently, Ai,i(0, x 2 ) = and A2,i(0, x\) — become the only possibility. 
Thus, dV\ is a subset of {x : X\ =0}. 

The same consideration applies to the pair A\ t2 {x) and A 2 ^ 2 (x) and shows that 
dV 2 is a subset of {x : x 2 = 0}. □ 

Corollary 3.10. Assume that A\, 2 — A 2 ^\ G nf . If WV is self-adjoint and the 
weight function W is consistent with T>, then 

Ai,i(x) = axx{c\ — xi), Ax i2 (x) = ax 2 (c 2 — xi), B 1 (x) = bx + d\) 

A 2 ,i{x) = axi(c3 — x 2 ), A 2 , 2 (x) — ax 2 (c 4 - x 2 ), B 2 (x) = bx + d 2 . 

Notice that (I3.14|l only imposes conditions on Aij on part of dV . If V is finite, 
there should also be conditions on the set dV \ {d\V U d 2 V). The equation (|3.13|) 
imposes such a condition, which is often fulfilled if the weight function W vanishes 
on the set dV \ {d x V U d 2 V). 

The property that d{V and d 2 V are subsets of {x : x 2 = 0} and {x : x\ = 0}, 
respectively, is fulfilled by the Type A and the Type B sets. In fact, they appear 
to be the only sets V that satisfy this property in connection with the difference 
equation. 

4. Difference equations and orthogonal polynomials 

4.1. Preliminary. The results in the previous section provides us with a method 
to identify difference equations that have orthogonal polynomials as solutions. The 
method consists of the following steps. 

Step 1. Among all admissible difference operators that satisfy H3.15|l . find T> that 
satisfies (|3.5f) . 
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Step 2. Using H3.6|) to identify the weight function that satisfies (|^ . 1|> . 

Step 3. Verifying that the condition i|3.13[l holds. 

We will follow the steps outlined above to do a case by case study. The classical 
discrete orthogonal polynomials of one variable will appear in the study, so we list 
them below (see QT] or [7]): 

Hahn polynomials u = Q n (x;a, f3, N): 

t(N +/3 + 1- «)AVu + (JV(a + 1) - {a + /3 + 2)t)Au = -n(n + a + (3 + l)u, 

where a, (3 > and TV e N. They are orthogonal with respect to the weight function 

mr;i on^o,!,...,^. 

Meixner polynomials u = M n (t; j3, c): 

tAVu + [c(t + (3) - t]Au = n(c - l)u 

where (3 > and < c < l. They are orthogonal with respect to the weight 
function ^c x on V = N . 

Krawtchouk polynomials u = K n (t;p,N): 

t(l - p)AVu + [p(N -t)- t(l - p)]Au = -nu 

where p > and N £ No- They are orthogonal with respect to the weight function 
O'tl-pr onV = {l,2,...,N}. 

Charlier polynomials u = C n (t;a): 

tAVu + (a - t)Au = -nu 

where a > 0. They are orthogonal with respect to the weight function a* jt\ on 
V = N . 

4.2. Difference equations with orthogonal polynomial solutions. We will 
now follow the steps outlined to identify the second order difference equations that 
have orthogonal polynomial as solutions. For this we assume that V is admissible 
and the coefficients Aij ofV satisfy H3.15J) . 

The equation (|3.5(l is a nonlinear equation, we solve it with the help of a computer 
algebra system. 

Case 1: a ^ 0. We can assume that a = — 1. Then Ai j and Bi must take the 
form a 

M,\{x) — x\(ci - xi), Ai^ 2 {x) = x 2 {c 2 - xi), Bi (x) = bx + di ; 
A 2 ,i(x) = xi(c 3 - x 2 ), A 2 ,2{x) = x 2 (c i - x 2 ), B 2 (x)=bx + d 2 . 

With the Ay and Bi given as above, solving the equation . 5|> (with the help of 
a computer algebra system) leads to essentially two solutions. 

Case 1 (i): One solution is 

d 2 = -{b + ci - c 2 )c 3 , di = c 2 (-b — c\ + c 2 ), c 4 = ci-c 2 + c 3 , 
so that, by solving (|3.6p. we get 

w( , _ r(-c 2 + xi)r(-c 3 + x 2 )r(ci + c 3 - xi - x 2 ) 

(X) ~ r(l + a:i)r(l + ar 2 )r(l + b + ci - c 2 - Xl - x 2 ) ' 
It is easy to verify that holds for V = {x : x\ > 0, x 2 > 0, x\ + x 2 < N}. 
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W a (x) 



This turns out to be exactly the Hahn polynomials of two variables discussed in 
0. Indeed, setting C2 = — (cti + 1), C3 = — (a 2 + 1), ci = N + a 2 + 0-3 + 2 and 
& = —3 — a 1 — (T2 — "3 shows that the difference equation is the one given below: 

Hahn Polynomials. 

(4.1) Xl (N -x 1 +(j 2 + 03 + 2)AiViu - x 2 (x! + o-i + l)AiV 2 « 

- xi(a;2 + ct 2 + l)A 2 Viu + x 2 (N - x 2 + a x + cr 3 + 2)A 2 V 2 u 

+ [(N - x 1 ){a l + 1) - x x (a 2 + <j 3 + 2)] A lU 

+ [(N-x 2 )(a 2 + 1) -x 2 {a x + 03 +2)] A 2 u = ~n(n + \a\ + 2)u, 

where \a\ = o\ + a 2 + 1T3. The weight function T4 7 is defined on V = {x : x\ > 
0,a;2 > 0,xi +x 2 < N} by 

^1 + <7l\ (X 2 + <72\ f N - X\ — £2 + 03^ 
Cl / V °"2 / V CT 3 

For each n, the equation has solutions 4>i >m {-\ cr), m = n — I and < I < n, given in 
terms of Hahn polynomials by 

4>i, m {x;a) = Qi{xf,a\,cr 2 + a 3 + 2m + l,N - l)(-N + x) m Q m (x 2 ; a 2 . a 3 , N - x\). 

Case 1 (ii): The other solution is 

di = c 2 (-b — cx + c 2 ), d 2 = -(l + b)c 3 , c 4 = l + c 3 , 
so that, by solving (|3.ti|l . we get 

, v = r(-c 2 +a:i)r(-l-6 + C2 + X2) 

1 ' r(i + Xl )r(i + x 2 ) 

The consistent condition (|3.13|) with g(x) = x$ shows that b = c 2 + c 3 — 1. Setting 
c 2 = —K\ — 1, c 3 — —k 2 — 1 and c\ = M + k 2 + 2 shows that the difference equation 
is of the same form as the equation (|4.1|l for the Hahn polynomials of two variables 
with (Ti = — «i, cr 2 — ~k 2 and = — 2 — k± — k 2 . Since M has to be a positive 
integer, this shows that K\ and k 2 need to be negative. However, in this case, 

w( s _ r(i + ki + xi)r(i + k 2 + x 2 ) 

[X '~ T(l + Xl )T(l + x 2 ) 

has singularities at xi = —K\ — 1 or x 2 — —k 2 — 1. Thus, this is not a proper 
solution. 

Remark. The computer algebra system that we used produced several solutions 
for 13. 5|) . upon close examination, however, we have (|4.1|l as essentially the only 
solution. The other solutions either turn out to be special cases of (|4.1|) or like the 
case (ii) given above. Even though the case (ii) does not lead to a proper solution, 
we have included the case to indicate how the method works. 

Case 2: a = 0. Assume that V is admissible and (|3.14|) holds. Then A^j and 
must take the form 

Ax t x{x) = a\xx, A\ t2 (x) — a 2 x 2l Bx(x) = bx + d\\ 
A 2 , 1 (x) = CL3X1, A 2 , 2 (x) = a i x 2 , B 2 (x)=bx + d 2 . 

With these Aij and Bi, solving the equation Q3.5J1 (with the help of a computer 
algebra system) shows that there are several solutions. 
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Case 2 (i): a 2 7^ and a 3 7^ 0. In this case we get 

di = a 2 d 2 /a 3 , ai — a 2 — b, 0,4 = 03 — b 
and, by solving (|3.6I) . 

(-a 2 r(-a 3 )y(a 2 + a 3 -b) N -*~y 



W{x) 



T(l - d 2 /a 3 - x - y)T(l + x)T(l + y) ' 

It is easy to verify that the condition (|3.13l) holds for V — {x : X\ > 0,x 2 > 
0, X\ + x 2 < N}. This turns out to be exactly the Krawtchcuk polynomials of two 
variables on V. Indeed, setting b = —1, a 2 = —p\, a 3 = —p 2 and d 2 = 0,3(1 — N) 
shows that the difference equation is the one given below: 

Krawtchcuk Polynomials. 

(4.2) (1 - pi)xi Ai Viit - pix 2 Ai V 2 u - p 2 xiA 2 Viu + (1 - p 2 )x 2 A 2 V 2 u 

+ [ Pl (N - xi) - (1 - pi)aui] Ai« + \p 2 (N - x 2 ) - (1 - p 2 )x 2 ] A 2 u 
= —nu. 

The weight function W is defined on V = {x : x\ > 0, x 2 > 0, X\ + x 2 < N} by 

For each n, the equation has solutions (f>i, m (-;p, N), m — n — I and < I < n, given 
in terms of Krawtchcuk polynomials by 

4>i, m {x; ,p,N) = Ki(x;pi,N - m)(—N + x) m K m (y;p 2 /(l - pi), N - x). 

These polynomials have appeared in J3| but the difference equation was not 
given there. 

Case 2 (ii): a 2 = and a 3 = 0. In this case we have several solutions, 
depending on the signs of various coefficients. These correspond to the product 
type orthogonal polynomials, which we list below. 

Product Meixner polynomials: 

(4.3) — i— -X1A1V1H+ 1 x 2 A 2 V 2 u 

Cl — 1 C2 — 1 

H [ci(xi + fix) - xi] Aim H 3_ [c 2 (x 2 + j3 2 ) - x 2 ] A 2 u = nu. 

Ci - 1 c 2 - 1 

The weight function W is defined on V = N§ by 

WbAx) = (M^ cT (M^ c? . 

Xl\ x 2 \ 

For each n, the equation has solutions </>/, m (-; 0, c), m = n — I and < I < n, given 
in terms of Meixner polynomials by 

<f> n (x;/3,c) = M l (x 1 ;l3 1 ,c 1 )M m (x2;(3 2 ,C2). 
Product Meixner-Krawtchcuk polynomials: 

(4.4) —^—rXx AiVitt - (1 - p)x 2 A 2 V 2 u 
c — 1 

H — [c(xi + (3) - xi] Aiu - [p(iV - a; 2 ) - (1 - p)x 2 ] A 2 u = mt. 
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The weight function W is defined on V = ■ i G No,0 < j < N} by 

For each n, the equation has solutions 4>i, m {'\ P, c,p, N) 1 m = n — I and < I < n, 
given in terms of Meixner and Krawtchcuk polynomials by 

<f) n (x;(3,c,p,N) = Mi(x 1 ;f3,c)K m (x 2 ; p, N) . 

Product Meixner- Charlier polynomials: 

(4.5) —^—x 1 A 1 \7 1 u - x 2 A 2 V 2 u 

c — 1 

+ 1 ^ [c(xi + /3) - xi] Aim - (a - x 2 )A 2 u = nu. 
The weight function VF is defined on V = Nq by 

si! x 2 ! 

For each n, the equation has solutions ^ jTO (-; /3, c, a), m = n — I and < / < n, 
given in terms of Meixner and Charlier polynomials by 

(p n (x;p,c,a) = Mi(xi;/3,c)C m (x 2 ;a). 

Product Krawtchcuk polynomials: 

(4.6) 

(1 - pi)xi Ai Vitt + (1 - p 2 )x 2 A 2 V 2 u 

+ [pi(iVi - xi) - (1 -_pi)xi] Aim + [p 2 (iV 2 - ar 2 ) - (1 - p 2 )x 2 ] A 2 u = -nu. 
The weight function W is defined on V = {(i, j) : < i < N u < j < N 2 } by 

x, /, , Nl - Xl fN 2 \ sN 2 -x 2 



w p (x) = ^ x i jpr(i-pir i - xi \ Kx *jp x 2 2 (i-P2 

For each n, the equation has solutions <^, m (-;p, N), m — n — l and < Z < n, given 
in terms of Krawtchcuk polynomials by 

cp n (x;p,N) = K t (x 1 ;p 1 , N 1 )K m (x 2 ;p 2 , N 2 ). 

Product Krawtcheuk- Charlier polynomials: 

(4.7) (1 -p)xiAiViu + x 2 A 2 V 2 u 

+ \p(N — x\) — (1 — p)x\] Aiu + (a — x 2 )A 2 u = —nu. 
The weight function W is defined on V = {(i, j) : < i < N, j G N } by 



2:2! 



For each n, the equation has solutions (f>^ m (-;a,p, N), m = n — l and < I < n, 
given in terms of Krawtchcuk and Charlier polynomials by 

<j) n (x;a,p,N) = Ki(x 1 ;p,N)C m (x2;a). 
Product Charlier polynomials: 

(4.8) X1A1V1U + x 2 A 2 V 2 u + (ai - x t )Aiu + (a 2 - x 2 )A 2 u = -nu. 
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The weight function W is defined on Nq by 

X\\ x 2 \ 

For each n, the equation has solutions ef>i m {-; a), m = n — I and < Z < n, given in 
terms of Charlier polynomials by 

(f> n (x;a) = Ci{x 1 ;a 1 )C m (x 2 ;a 2 ). 

Case 2 (iii), other possibilities. There is no solution if one of the a 2 or 03 is 

zero and the other one is not. However, if we set 04 = 0, then solving the equation 
(pH)|) leads to 

di = — (ai + b)d 2 /b, a 2 = ai + b, 03 = b. 

Solving the equation (|3.6[l leads to the weight function 

W(x) = b y ( ai + b) N - y T(d 2 /b + x + y)/(xly\). 

Since W does not have enough decay, the set V needs to be bounded to keep the 
linear functional C defined by W well-defined on polynomials. However, we do not 
see a way to choose V such that W vanishes on the boundary dV \ {d\V H d\V). 
In other words, we do not see a way to choose V such that the condition l|3.13[) can 
be satisfied. 

4.3. Final Remarks. Under the assumption that A i: j are of the form (|3.15|) . we 
found eight difference equations that have orthogonal polynomials as solutions. 
These are the Hahn and the Krawtchcuk polynomials on Type A set V, and product 
type polynomials on Type B set V: product Meixner, product Krawtchcuk, and 
product Charlier polynomials, as well as product Meixner-Krawtchcuk, Meixner- 
Charlier, and Krawtchcuk-Chalier polynomials. With the restrictions that we put 
on the difference equations (that is, the conditions in Corollary 3.10), these appear 
to be the only possible solutions. 

Let us mention that another attempt of characterizing orthogonal polynomials 
of two variables satisfying second order difference equations, based on a matrix 
approach, is currently being undertaken by a group in University of Granada f|l(J|). 

As pointed out in the introduction, if u satisfies a difference equation then the 
function Uh(x) := u(h±xi, h 2 x 2 ) satisfies a similar difference equation which be- 
comes in limit, as h — > 0, the second order partial differential equation. The second 
order partial differential equations that have orthogonal polynomials as solutions 
are characterized in jHj. It is known that there are nine such equations, whose 
solutions correspond to product type orthogonal polynomials of two variables, as 
well as two other type, one is the Jacobi type orthogonal polynomials on the tri- 
angle T = {(x, y) : x > 0, y > 0, < x + y < 1} and the other is the orthogonal 
polynomials on the disk B = {(x,y) : x 2 +y 2 < 1}. The discrete Hahn polynomials 
of two variables on Type A set appear to become the Jacobi type polynomials on 
the triangle. However, as far as we can tell, there is no discrete analog for the 
orthogonal polynomials on the disk. 

A natural question is if there are other difference equations of the form (|1.3|) 
that have orthogonal polynomials as solutions. Notice that our results are obtained 
under various assumptions on T>. For example, we assume that the quadratic parts 
of Ai >2 and A 2t i are equal, that is, Ax >2 — A 2 ^\ e H\ . Furthermore, we establish 
the orthogonality by requiring that V is self-adjoint and W is consistent with V. 
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These assumptions are sufficient but may not be necessary. Still, assuming that 
Ai,2 — Aq.,1 G n^, it seems to us that likely no other nontrivial difference equations 
can have orthogonal polynomials as solutions. 
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